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Abstract 

Magnetization processes and phase transitions in a geometrically frustrated 
triangular lattice Ising antiferromagnet in the presence of an external magnetic 
field and a random site dilution are studied by the use of an effective-field theory 
with correlations. We find that the interplay between the applied field and the 
frustration-relieving dilution results in peculiar phase diagrams in the temperature- 
field-dilution parameter space. 
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1 Introduction 



The Ising model with nearest-neighbor antiferromagnetic interactions on a tri- 
angular lattice is fully frustrated due to local geometric constraints that prevent 
simultaneous minimization of all the pairwise interactions. It has been solved ex- 
actly and found to display no ordering at any finite temperatures [TP] . Only at T = 
K the model displays a critical phase with algebraically decaying correlations [5]. 
The ground state is highly degenerated and thus even a small perturbation can cause 
that the system will choose an ordered state. Such a perturbation can be an external 
magnetic field, which by several studies has been shown to produce a line of second- 
order phase transitions for a certain range of the field values |4H6] . At the transition 
the system passes from the ferrimagnetic phase with two sublattices aligned parallel 
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and one antiparallel to the field at lower temperatures (itt) to the paramagnetic 
phase in which all spins are aligned parallel to the field at higher temperatures (ttt) • 
As a result of the frustration, at low temperatures the magnetization versus field 
curves display a broad plateau with the value of m = 1/3. Similar stepwise plateaus 
with different magnetization values have also been observed in some other lattices, 
such as the recursive lattice [7], the zig-zag ladder |8l|9] and even on the cluster 
level p^fTT]. However, besides the magnetic field, the degeneracy of the system can 
also be removed by considering the interactions of more distant neighbors p!2] - [T5] . 

Diluting the system with nonmagnetic impurities in zero field locally relieves 
frustration and can lead to a spin-glass order [ISIIIZ], although this scenario re- 
mains controversial for the current 2D system [18]. Spin-glass ordering has also been 
observed in another highly frustrated Ising antiferromagnet on a fee lattice for a cer- 
tain range of the magnetic concentration [T9]. In Ref. |20j the effects of doping the 
triangular Ising antiferromagnet with slow moving holes has been studied. Besides 
the site-dilution problem, also a bond-dilution problem and a random-bond 
Ising model [22] on a triangular lattice have been investigated. In contrast to the 
uniform dilution, if only one sublattice of the triangular lattice is diluted then a 
long-range order can develop in the remaining two sublattices already at relatively 
low concentrations of the impurities [23] • Very recently, Yao [21] has studied how 
dilution can modulate the frustration effect in the system in a field. For example, 
he found that even a small dilution causes that the broad frustration-induced 1/3 
magnetization plateau splits into a stepwise curve. The importance of such studies is 
amplified by the fact that there are some corresponding real magnetic compounds, 
such as Ca3Co206, that have been experimentally found to display some unconven- 
tional frustration-induced features [22] and the current model has been successfully 
employed in their explanation [26] [27]. 

Motivated by the previous interesting findings, the objective of the present work 
is to systematically study effects of the uniform site dilution on the phase transitions 
in the triangular lattice Ising antiferromagnet in a field. 



2 Model and effective-field approach 

The model Hamiltonian is given by 

H=-JY: UjSiS,-hJ2^^S^1 (1) 

<i,j> i 

where Si = ±1, are the Ising spin variables, h is the external magnetic field, J < is 
the exchange interaction constant, and < i,j > is the sum extending over all nearest 
neighbor (NN) pairs, are quenched, uncorrelated random variables chosen to be 
equal to 1 with probability p, when the site i is occupied by a magnetic atom and 
0, with probability 1 — p otherwise. Then the probability distribution is given by 
P{^i) = pS{^i — 1) + (1 — and p represents the mean concentration of magnetic 
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sites. 

We employ an effective- field theory (EFT) with correlations (see e.g., [28]), based 
on a single-site cluster approximation with the attention focused on a cluster com- 
prising just a single spin, labeled i, and the NN spins with which it directly interacts. 
Following Ref. [29] we get the exact relation 

US,) = Utanh[^{jj2^A + h)]), (2) 
i=i 

where z is the number of NNs of the site i (i.e., the coordination number), /3 = 
l/ksT, and (. . .) denotes a thermal average for a fixed spatial configuration of the 
occupied sites. Applying the differential operator technique [3D] to the identity ([2]) 
and using the exact relations 

exp(A^j) = exp(A) + 1 - ^j, (3) 

exp(/iS'j) = cosh(/i) + Sj sinh(/i), (4) 

one obtains 

USi) = d nfe- cosh(^JL>) + ^jSj sinh(/3JD) + 1 - ^j]) tanh{x + /3/i)U=o, (5) 

where D = d/dx is the differential operator. In order to carry out the configurational 
averaging over the occupational numbers ^i, let us assume that NNs of the site i are 
completely independent of each other by taking an approximation 

{S,S,...Si)^{S,){S,)...{Si). (6) 

In spite of this simplification, we note that this approximation is quite superior to 
the standard mean-field theory, since here, by using Van der Waerden identity dl]), 
the relations like Sj=l are exactly taken into account. 

However, unlike on some other lattices, on a triangular lattice NNs of the central 
spin i include pairs of spins that are also mutual NNs and, therefore, their decoupling 
by the approximation ([6|) could result in rather high inaccuracies. Even more impor- 
tantly, such a straightforward application of EFT would lead to a complete lost of 
the frustration and therefore inevitably incorrect results. In order to include all the 
NN interactions and the effect of the geometrical frustration, we chose to partition 
the lattice into three interpenetrating sublattices A, B and C in such a way that 
spins on one sublattice only interact with spins from the other two sublattices (see 
Fig. [1]). Then all the NN interactions are accounted for and the frustration arises 
from the effort to simultaneously satisfy all the mutual intersublattice interactions, 
which are antiferromagnetic and isotropic i.e.. Jab = Jac = Jbc = J < 0. As an 
example, let us consider the situation in a selected triangular plaquette i in Fig. [1] If 
the spin SiA is in the state +1 then the energy is minimized if all its NNs, including 
SiB and Sic, are in the state —1. However, for the spin SiB in the state —1 it would 
be energetically favorable if all its NNs, including S^a and Sic, were in the state 
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which creates frustration for the spin Sic- Given the lattice partition, the expression 
([5]) can be written for individual sublattices in the form 

US^) =d n cosh(/3 JD) + ^,S, smh{f3JD) + 1 - Q 

X n [6 coshiPJD) + ^kSk sinh(/3JD) + 1 - a]) tanh(x + /3/i)U.=o, (7) 
k=i ' 

where zi, Z2 are the numbers of NNs of the spin Si on a given sub lattice that belong 
to the remaining two sublattices. Then, by performing configurational averaging 
for all three sublattices, the respective sublattice magnetizations mx = {C,i{Six))c, 
X=A,B and C, can be calculated from the set of coupled equations 

niA = p{a + bniB)^ (a + hmcf tanh (x + (3h) \x=o^ 

m-Q = p{a + hmxf (a + bmcf tanh (x + (ih) \x=o, (8) 

niQ = p (a + hrax)^ (a + hm-of tanh [x + I3h) \x=q-, 

where a = 1— p+pcosh {PJD), b = sinh {(5JD) andp = {C,i)c is the same for all three 
sublattices. The explicit form of Eqs. ([8]) can be calculated by using the mathematical 
relation exp{aD)f{x) = f{x + a). Then we can define the total magnetization per 
site m = (mA + rriB + mc)/3 and the order parameter o = [max{m\,mB,Tnc) — 
min(mA, "niB, '^c)]/2. The latter serves to localize phase boundaries between the 
ferrimagnetic (J, ft) the paramagnetic (ttt) phases. 

[Fig. 1 about here.] 



3 Results and discussion 

In zero field, Eqs. ([8]) have only trivial solution at all temperatures, which means 
that our effective-field approach reproduces the exact result of no long-range order 
down to T = K [Tl[2] . In Fig. [2] we plot the field dependence of the magnetization 
at low temperatures for the pure case (p = 1) as well as diluted cases {p < 1) 
with different degrees of dilution. For the pure case, the curve displays the typical 
frustration- induced broad 1/3 plateau, in accordance with some previous theoretical 
[2HEI] and experimental [22] observations. When the system is diluted, we can 
observe formation of multiple steps at about integer values of the field h/\J\ = 
1, 2, 3, 4, 5, which get more pronounced as the dilution increases. This behavior is in 
agreement with the Monte Carlo results [21] and it could be attributed to a gradual 
flipping of some dilution-relieved spins to the field direction when the Zeeman term 
contribution overcomes their exchange energy. 

[Fig. 2 about here.] 
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[Fig. 3 about here.] 



In order to establish the phase diagrams in the temperature-field-dilution param- 
eter space, we monitor the value of the order parameter o. Namely, the transition 
temperature is set to the value at which the order parameter vanishes. For p = 1, 
the phase boundary (curve (a) in Fig. [3]) agrees quite well with those obtained in 
the previous Monte Carlo mean-field, Monte Carlo and renormalization-group stud- 
ies [1H6]. The quenched dilution in the present system relieves frustration, however, 
only locally and thus does not support long-range ordering, in contrast to the se- 
lective dilution case Indeed, as evidenced in Fig. E] (see curves (b)-(i)), the 
increasing dilution lowers the transition temperatures. Furthermore, it shrinks the 
region of the field values at which the transition can occur. As already signaled by 
the magnetization behavior shown in Fig. [21 the shrinking occurs discontinuously 
and asymmetrically. More specifically, the initial interval of the values at which the 
transitions can occur h/\J\ G (0, 6) for p = 1 shrinks discontinuously to (0, 5) at the 
concentration pi = 0.96, then to (1,5) at p2 = 0.941, etc. until pq = pc = 0.802, 
below which no long-range order can survive. Such behavior is illustrated in Fig. HI 
where we plot the sublattice magnetizations, the total magnetization, and the order 
parameter as functions of the applied field for two values of p = 0.95 > p2 and 
p = 0.94 < p2, at kBT/\J\ = 0.1. 

[Fig. 4 about here.] 

Owing to the fact that the phase boundaries always terminate at integer val- 
ues of the field as the ground state is approached, for some values of the field and 
concentration we can observe a reentrant phenomenon. For example, it occurs for 
p = 0.95 (curve (b) in Fig. |3]) at /;,/| J| = 5.05, as shown in Fig. |5l The system passes 
from the paramagnetic phase at low temperatures to the ferrimagnetic one at higher 
temperatures and back as temperature is further increased. Such a reentrant behav- 
ior can be observed in a low- field region (just below h/\J\ = 1) at concentrations 
p G (0.932, 0.941) and in a high-field region (just above h/\J\ = 5) at concentrations 
p G (0.934,0.960). Thus for p G (0.934,0.941) the reentrant phenomena appear at 
both ends of the boundary curve. 

[Fig. 5 about here.] 

In Fig. [3] we can also observe that as dilution is increased the boundary curves 
start to develop a valley. First it appears forp = 0.904 at the fields just above h/\J\ = 
4 (see curve (d)) and gets deeper as dilution increases (curve (e)) until it disappears 
at p = 0.876 and the phase transition field interval shrinks from (1, 5) to (1, 4) (curve 
(f)). At the same time, another valley emerges at the fields just below h/\J\ = 2 and 
persists from p = 0.883 down to p = 0.865, below which the phase transition field 
interval shrinks again from (1,4) to (2,4). Eventually, below p = 0.848 yet another 
valley appears close to h/\J\ = 3 (curves (g),(h)) until it disappears at p = 0.833 
and the phase boundary becomes limited to the fields h/\J\ G (2,3) (curve (i)). 
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The presence of multiple valleys in the same {kBTc/\J\,h/\J\) phase boundary for 
some values of p enables us to observe reentrant phenomena by varying the field at 
a fixed temperature. In particular, a typical pure system behavior of passing from 
the paramagnetic phase at low fields to the ferrimagnetic one at higher fields and 
back can change to such that, depending on the dilution and temperature values, 
the system can repeat this cycle twice or even three times, as shown in Fig. [61 
The latter situation can be observed for p G (0.877,0.879). As a result, the system 
displays six percolation thresholds pj, i = 1, . . . , 6, corresponding to different field 
intervals, as listed in Tabled] The phase diagram in the p — h/\J\ plane close to the 
ground state (at ksT /\J\ = 10~^) is shown in Fig. [71 The lines exceeding the edges 
at h/\J\ = 2,3,4 signify persistence of the reentrant phenomena described above. 

[Table 1 about here.] 

[Fig. 6 about here.] 

[Fig. 7 about here.] 



4 Conclusions 



We have applied an effective-field theory with correlations to study phase tran- 
sitions in a geometrically frustrated site-diluted triangular lattice Ising antiferro- 
magnet in a field. For the pure system our approach reproduced the exact solu- 
tion of no long-range order in zero field [Tl|2] and in a finite field yielded a good 
agreement with the results obtained by more accurate approaches, such as Monte 
Carlo simulations [1] and renormalization-group theory [S]. To our best knowledge, 
we are not aware of any attempts to establish phase diagrams of the present di- 
luted system. However, our results presented in Fig. [2 nicely corroborate the recent 
Monte Carlo density of states investigations of magnetization processes in such a 
system [23], namely the stepwise splitting of the low-temperature broad 1/3 mag- 
netization plateau at integer values of the field as a result of dilution. Although 
in [23] no order parameter was evaluated and therefore it is not clear which of the 
jumps in the magnetization are associated with the phase transitions, it is evident 
that the low-temperature transitions in that study should occur at some integer 
values between < h/\J\ < 6, in accordance with our results. Since the peculiar 
critical behavior with multiple reentrants observed in our study is a consequence 
of such discontinuities, we speculate that they are not merely artifacts of the used 
approximation and might be also reproduced in Monte Carlo simulations. Notwith- 
standing, further theoretical and experimental evidence is desirable to confirm our 
predictions. Suitable candidates for the experimental studies could be compounds, 
such as CagCoaOe or CS3C0X3, (X=C1 or Br) [25]. 
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Fig. 1. Triangular lattice partition into three sublattices A,B and C. 
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Table 1 

Percolation thresholds corresponding to different field h/\J\ intervals. 

h/\J\ (0,1) (1,2) (2^3) (3,4) (4,5) (5,6) 

Pi p2 = 0.941 p4 = 0.865 p6 = Pc = 0.802 p^ = 0.833 pa = 0.876 pi = 0.96 
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